Dynamic conductivity in graphene beyond linear response 
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The independence of the dynamic conductivity of intrinsic graphene of frequency takes its origin 
in the compensation of the vanishing density of states by the diverging matrix element of the 
corresponding interband transition. The applicability of the linear response approach, however, 
breaks down when this matrix element becomes comparable with the inverse electron lifetime. We 
show that the physics of the ac conductivity in this regime is determined by Rabi oscillations and 
obtain it beyond the first order perturbation theory. Under strong applied electric fields the induced 
current eventually saturates at a value determined by the frequency and the lifetime. We also 
calculate the electromagnetic response of a graphene sheet and find that the optical transparency is 
increased by the non-linear effects and make experimental predictions. 

PACS numbers: 73.23.-b, 72.25.+i 



Introduction. Low-temperature dynamic conductivity 
(fcflT <C TwS) of graphene was predicted to be of the or- 
der of the conductance quantum and independent of fre- 
quency, (To = c 2 jAJi 1], as long as the frequency ui is small 
compared with the electron bandwidth [2fl. This predic- 
tion has been verified experimentally 

®MM- Theoret- 
ical attempts to understand the role played by electron- 
electron interaction began even earlier [6(] . Following the 
conclusion that the interactions do not modify the con- 
ductivity in the limit u> — > 0, [1| it was further pre- 
dicted that even at non-vanishing frequencies the inter- 
action corrections are extremely small 0, and should 
hardly be noticeable at all, the claim consistent with re- 
cent precision measurements of the optical transparency 
in the visible frequency range 

In this Letter we address a different aspect of the phe- 
nomenon of the minimal conductivity which has not been 
discussed before. We begin by pointing out that the 
minimal conductivity is essentially a linear response con- 
cept, which needs to be reevaluated in relation to Dirac 
fermions. For comparison, in a conventional metal the 
linear response formalism is reliable as long as the energy 
acquired in the external electric field over one period of 
its oscillation, eEv/u, (or the mean free time, if disorder 
scattering can not be ignored) must be negligible com- 
pared with the Fermi energy Ef that is typically very 
large. 

In contrast, in clean intrinsic graphene Fermi energy 
vanishes, Ep — 0, and the above reasoning fails. To 
gain a better understanding of the essential physics, let 
us take a closer look at the Golden rule derivation of 
the ac conductivity do- Fig- 1 illustrates a vertical in- 
terband transition between the two states in the filled 
lower cone and the empty upper cone, which is responsi- 
ble for the energy dissipation in electric field. The matrix 
element of this transition is of the order of the energy 
gained in external field, H12 ~ eEv/uj. The number of 
available electrons in the lower cone with energy Klo/2 
is determined by the density of states which is linear in 
frequency, D(u>) ~ u>/v 2 h. Golden rule expression for 



the energy absorbtion rate in these transitions is sim- 
ply, £ ~ uj\H 12 \ 2 D(lu) ~ (e 2 /h)E 2 . Recalling the usual 
Joule heat expression we arrive at the electric conduc- 
tivity independent of frequency and of the order of the 
conductance quantum. 
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FIG. 1: Interband transitions induced by external electric 
field of frequency uj. Dashed lines stand for a single transi- 
tion (short lifetime 1/r) or multiple Rabi oscillations (long 
lifetimes 1/r). Dotted lines illustrate energy dissipation from 
relaxation of electrons in the upper Dirac cone and holes in 
the lower cone. 

It is well-known, however, that the first order pertur- 
bation theory fails close to a resonance, i.e. when the 
matrix element becomes of the order of the detuning en- 
ergy, which in our problem is Huj — 2vp. Since available 
momenta are continuous there are always states with the 
vanishing detuning, in which case the inverse lifetime T, 
calculated for a state with the corresponding momentum, 
must be used as a cut-off. This condition imposes the re- 
striction upon the strength of electric field, eE <§; TiwY/v. 
As frequency decreases, u> — > 0, the amplitude of applied 
field must vanish significantly faster, at least as fast as 
~ lo 2 12]. Here we concentrate on the response at the 
same frequency as the frequency of external field, disre- 
garding possible effects of frequency multiplication. From 
the above arguments one can expect that the amplitude 
of induced current will be, 



j =o- Q E[l + C 



e 2 E 2 v 2 
h 2 uu 2 T 2 
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and singular at low frequencies. In the present Letter we 
identify the leading terms in the series |T]) with the Rabi 
oscillations between states in the lower and upper cones 
with the same momentum. 

The above analysis indicates that the lifetime T plays 
a central role in the calculation of ac conductivity. When 
r is small, or more precisely, when the dimensionless pa- 
rameter 



c 



evE 
TiujT 



(2) 



is large, Rabi oscillations between states in the two cones 
persist for a long time and make energy dissipation less 
effective thus decreasing the conductivity. In the oppo- 
site limit, £ <C 1, fast relaxation ensures that the states 
in the upper cone are always empty and that no Rabi os- 
cillations can occur, thus making first-order perturbation 
theory applicable. We find that a homogeneous external 
electric field applied along the x-axis within the plane of 
graphene 



E(t) = ±Ecos ut. 



(3) 



induces the electric current which for arbitrary values of 
( has the form (l3| . 



2<tqE cos tot 



(4) 



We observe that the minimal conductivity is recovered 
when C — > 0, however for strong electric fields (or small 
frequencies), when ( ^> 1, the current saturates at the 
value 



3r, 



\e\ujT 
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exactly the same momentum. Denoting by a p (t) the am- 
plitude to occupy state in the lower cone and by b p (t) the 
same for the upper cone we write the total wave function 
as. 
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(8) 



where Xp is the angle between the electron momentum p 
and the x-axis. Spinor Schrodinger equation id^ p /dt = 
Hipp reduces to a pair of equations (we assume for brevity 
that h — 1), 



ib p = 



evEo 



evEg 



sinujt (ib p sinXpf 



-2ivpt 



a P cos x P ), 



sin tot (bp cos Xp — *flp sin Xpe 
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(9) 



In solving Eqs. © we follow the established route 
by keeping only the low-frequency resonant terms, ~ 
e ±i{uj-2vp)t ^ £ ne approximation known in quantum op- 
tics as the rotating wave approximation. An additional 
modification of the equations of motion is needed, how- 
ever. In the present form Eqs. ([9]) lack relaxation and 
describe undamped Rabi oscillations leading to the van- 
ishing energy dissipation. To include energy relaxation 
we take into account the finite lifetime l/r p , which due 
to electron-hole symmetry is assumed to be the same for 
electrons in the upper cone and holes in the lower cone: 



-^ pe ^- 2 ^) t 6 p + ^(l-a p ), 



-idpe-^-^ap - 



(10) 



This result can be qualitatively visualized as one electron 
being transferred per its lifetime 1/r across each strip of 
graphene that has the width equal to the electron wave- 
length v/uj corresponding to energy Tilu. 

Derivation. Homogeneous external electric field is 
most conveniently incorporated in the electron dynam- 
ics with the help of the longitudinal gauge and vector 
potential 



A(t) — ~x—E smut. 

LU 



(6) 



The Hamiltonian of graphene is the well-known Dirac 
Hamiltonian in the pseudospin (sublattice) space 



H = vcr ■ p H o x E smut. 

LU 



(7) 



where we introduced the Rabi frequency fi p = 
^^sinxjj- While the relaxation is included into 
Eqs. (|10p in a phcnomcnological way, it has an advantage 
of allowing analytical solution. After straightforward cal- 
culations, 



a p (t) 



b P (t) = 



rg/4-»r p A p /2 
fi| + T2/4 - iT p A p /2 
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fi p + I*/4 - zT p A p /2 
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The crucial observation that allows the solution of the 
problem beyond first order perturbation theory is that 
the momentum p is an integral of motion. The exter- 
nal electric field therefore mixes only the two states with 



where A p = lo — 2vp denotes the detuning between the 
frequency of electric field and the electron-hole excitation 
energy. The terms containing constants C+ and C_ de- 
pend on the initial conditions but decay with time and 
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thus are unimportant for the properties of a "stationary" 
state. 

Expectation value of electric current carried by a pair 
with momentum p is given by, jp(t) = ev^(t)a x ^ p (t). 
Note that the vector potential term in Eq. (0 does not 
modify the expression for current since the former is in- 
dependent of momentum. We obtain, 

e 2 v 2 Tlcosujt + 2TlA p smujt „ / x 

m= E 8uj ,l „,,2 p i.,,. sin x P - (ii) 



The total current is found upon summation over all 
momenta, two spins directions and both Dirac points, 
j = 4^ p j p . For the first contribution, ~ cos cut, which 
is responsible for energy dissipation, it gives, 



jx {t) = cos Lot 



'-E 



2-7T OO 



2tt 2 lu 



o o 



T p pdp sin 2 xdx 



(12) 



where we introduced the notation 



e 2 vE 2 sin \ 

4T p LJ 2 ' 



Provided that f3(p,x) the integrand is a sharply 

peaked function of p — lo/2v and the integration over dp 
can be carried out easily. Two conditions must be sat- 
isfied in order to ensure this. The first one is essentially 
the condition that electrons are well-defined excitations, 
r«w. The second condition imposes the restrictions on 
the upper value of electric field. In terms of the dimen- 
sionless parameter defined earlier, Eq. ((2]), this condition 
implies that £ <C lo/T. 

Integrating the sharply peaked Lorentian over dp and 
taking the subsequent integration over the angle we ob- 
tain, 



jx (t) = cos 




(13) 



which is equivalent to Eq. ([4]). The lifetime in Eqs. ([2][4]) 
and (fl~5]) is taken at the value of momentum correspond- 
ing to frequency lj: T = T^/ 2v . 

The second term in the numerator of Eq. pip describes 
the out-of-phase non-dissipative part of the current. This 
term is odd in A p = lj — 2vp and thus vanishes after 
integration over the vicinity of resonant momenta. Note 
that the contribution of non-resonant momenta to the 
~ svciLot term though small is not necessarily zero, but is 
beyond the reach of the approximations employed in our 
analysis. 

Scattering mechanisms. The dependence of the inverse 
lifetime T(u>) on frequency determines the value of the 
saturation current (O and depends on the importance of 
various scattering mechanisms that could be present in 
graphene. 



i) In case of short-range impurities (vacancies, non- 
charged substitutions) the electron lifetime in the first 
Born approximation is simply proportional to the density 
of states, T(lj) cx D(u>) cx u> and the saturation current 
grows as the second power of the frequency of electric 
field. 

ii) When charged impurities dominate, the amplitude 
of the bare Coulomb impurity potential is inversely pro- 
portional to the electron momentum, 2ire 2 /q cx 1/lu. 
Since the screening in intrinsic graphene does not mod- 
ify this dependence, the lifetime T(lo) cx D(lo)/lo 2 cx 1/lj 



is inversely proportional to the frequency [14{ . Thus, 
the value of the saturation current is independent of fre- 
quency. In both cases of neutral and charged impurities 
jmax is inversely proportional to the concentration of im- 
purities. 

iii) Electron-phonon scattering at frequencies above 
the Debye frequency is quasi-elastic and follows the same 
linear dependence originating from the density of states 
as is the case for short-range disorder, T cx 1/lo [lj|. 

Implications for optical experiments. Of particular sig- 
nificance for ongoing experiments [H, 0, [H, El is the cal- 
culation of the response of a graphene layer to an incident 
electromagnetic wave in the nonlinear regime. Below we 
assume graphene suspended in vacuum (or air) and de- 
termine its optical transparency. 
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FIG. 2: Dependence of the transmission coefficient T of sus- 
pended graphene on the effective intensity of the incident ra- 
diation, I = (evEo/huiF) 2 = At low intensities it starts 
from the value of 97.7% corresponding to the first order per- 
turbation theory but increases sharply with the intensity I, 
reaching 99% for (o ~ 3. 



The incident radiation is assumed to be propagating 
along the negative direction of the z-axis and normal to 
the plane of graphene (xy-plane) with the linear polar- 
ization directed along the x-axis. The electric field con- 
sists of the incident Eq and reflected En waves above the 
graphene flake, and of the transmitted wave Et below it 



E : 



E cos (ujt +^f)+E R cos (ujt - *f 



■), *> 0.(14) 
z < 0. v ; 
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Maxwell's equations with the two-dimensional current 
provided by graphene sheet at z = taken into account 
has the form, 



d 2 E 1 d 2 E 2na a uj E T smwt 



dz 



dt 2 



5(z), (15) 



where (t = corresponds to the field acting in the 

graphene flake, Et- The first boundary condition re- 
quires the continuity of the electric field, E$ + En — 
Et- The second condition is obtained from integrating 
Eq. (fTS"|) across the plane of graphene. 

Solving for Er and Et we obtain the following al- 
gebraic equation for the transmission coefficient T = 
E 2 /E 2 



0! 



na 



1 + VI + TI y/f' 



(16) 



where a = e 2 /Tic stands for the fine structure constant, 
and the dimcnsionlcss intensity of the incident wave is 
introduced, 7 = (evE /hujT) 2 . Fig. [2] shows the depen- 
dence of the transmission coefficient on the magnitude of 
the incident wave. At low intensities, 7 — > the trans- 
mission approaches the linear response value T = 97.7%, 
reported in the experiment ll|. As the intensity grows 



the transmission initially increases steeply, but eventu- 
ally tends much more slowly towards T = 1 according 
to, T pa 1 - 27ra/V7. 

The modern experimental capabilities are well within 
the necessary means to resolve the nonlinear transmis- 
sion, as predicted by Fig. [5] In particular, the preci- 
sion reported in Ref. [ll[ was 0.1%. It is interesting to 
note the upward trend in T(A) in the data presented 
in Ref. [llj for different wavelengths A. It is qualita- 



tively reminiscent of the T(7) dependence predicted in 
the present Letter, if one assumes short-range impurities 
and/or phonons as the primary scattering mechanisms. 
Understandably, systematic and more detailed measure- 
ments of the intensity dependence are necessary before 
any comparisons could be made with the theory. 

Summary. The conductivity of intrinsic graphene, 
e 2 /Ah, has been established to remains constant even 
when the frequency tends to zero despite the vanishing 
of the density of available states, D(u>) oc u — > 0. This 
celebrated result originates in the fact that the matrix 
element, evE/uj, for the corresponding interband transi- 
tion between a filled state in the lower Dirac cone and an 
empty state in the upper cone diverges at low frequencies. 
This divergence implies the failure of the first order per- 
turbation theory and, therefore inapplicability of the min- 
imal conductivity at low frequencies for a fixed strength of 
applied electric field. The Golden rule approach is valid 
only if the magnitude of electric field vanishes at least as 
uj 2 , so that (< 1. In this paper we extended the theory 
of electrical conductivity towards finite electric fields for 
arbitrary values of C, as long as C *C lu/T(uj). The key 



observation that allows a solution beyond perturbation 
theory is the fact that momenta are conserved since the 
important transitions are direct, see Fig. |T]). The prob- 
lem is therefore reduced to strong coupling within sim- 
ple (and independent of each other) two-level systems. 
The electric conduction is largely determined by those 
states that are in the vicinity of the resonant transition, 
|cl>— 2vp\ ~ r. When the matrix clement for the transition 
exceeds the inverse lifetime 1 /T multiple Rabi oscillations 
develop that greatly reduce the probability of energy dis- 
sipation (Joule heat) and, hence, conductivity. As shown 
above this must result in the experimentally observable 
enhancement of the transparency of graphene. 

The existing experiments have been performed in the 
infrared [1, 0, [B| and visible light 



ll| domains. Previ- 



ously, however, nonlinear effects have not been systemati- 
cally searched for. We expect that the stronger nonlinear 
signatures could be revealed for clean samples, strong in- 
tensities of incident radiation and lower frequencies (pos- 
sibly infrared). 
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ported by DOE, Office of Basic Energy Sciences, Award 
No. DE-FG02-06ER46313. 



[1] 
[2] 



[3] 

[4] 

[5] 

[6] 
[7] 

[8] 

[9] 
[10] 

[11] 

[12] 
[13] 



A.W.W. Ludwig, M.P.A. Fisher, R. Shankar, and G. 
Grinstein, Phys. Rev. B 50, 7526 (1994). 
In the present paper we consider clean graphene and as- 
sume that frequency is larger than the broadening of elec- 
tron states, lo r. In the opposite limit, w <C L, the 
linear conductivity has to be determined by means of a 
self-consistent theory, see e.g. S. Ryu, C. Mudry, A. Fu- 
rusaki, and A.W.W. Ludwig, Phys. Rev. B 75, 205344 
(2007), and is numerically different, a — 2e 2 /ty 2 Ti. We 
leave the nonlinear effects in a for future study. 
Z.Q. Li, E.A. Henriksen, Z. Jiang, Z. Hao, M.C. Martin, 
P. Kim, H.L. Stormer, and D.N. Basov, Nature Physics 
4, 532 (2008). 

A.B. Kuzmenko, E. van Heumen, F. Carbone, and D. 
van der Marel, Phys. Rev. Lett. 100, 117401 (2008). 
K.F. Mak, M.Y. Sfeir, Y. Wu, C.H. Lui, J.A. Misewich, 
and T.F. Heinz, Phys. Rev. Lett. 101, 196405 (2008). 
E.G. Mishchenko, Phys. Rev. Lett. 98, 216801 (2007). 

D. E. Sheehy and J. Schmalian, Phys. Rev. Lett. 99, 
226803 (2007). 

I.F. Herbut, V. Juricic, and O. Vafek, Phys. Rev. Lett. 
100, 046403 (2008). 

E. G. Mishchenko, Europhys. L ett. 83, 17005 (2 008). 
D.E. Sheehy and J. Schmalian. larXiv:0906.5T64l 

R.R. Nair, P. Blake, A.N. Grigorenko, K.S. Novoselov, 
T.J. Booth, T. Stauber, N.M.R. Peres, and A.K. Geim, 
Science 320, 1308 (2008). 

For the quasiparticle picture to be meaningful it is in any 
case necessary to have T <C lo. 

In this paper we concentrate mostly on the part of 
the current that is in phase with external electric field, 
j(t) — j cos tot rather than the out-of-phase contribution, 



5 



~ sinujt, which does not lead to energy dissipation. See [15] C-H. Park, F. Giustino, M.L. Cohen, and S.G. Louie, 
also explanation below Eq. ([13}. Phys. Rev. Lett. 99, 086804 (2007). 

[14] K. Nomura and A. H. MacDonald, Phys. Rev. Lett. 98, 
076602 (2007). 



